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ABSTRACT 



K*" I The MACHO project has been monitoring about ten million stars in the 

I Large Magellanic Cloud in the search for gravitational microlensing events 

^ _ caused by massive compact halo objects (Machos) in the halo of the Milky 

O ! Way. In our standard analysis, we have searched this data set for well sampled, 



long duration microlensing lightcurves, detected several microlensing events 
consistent with Machos in the 0.1 Mq ^ m < 1.0 M© mass range, and set 



Qh| limits on the abundance of objects with masses 10 Mq < m < 10 ^M©. In 

O 
H 
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this paper, we present a different type of analysis involving the search for 
very short time scale brightenings of stars which is used to set strong limits 
on the abundance of lower mass Machos. Our analysis of the first two years 
of data toward the LMC indicates that Machos with masses in the range 
2.5 X 10"'^ Mq < m < 5.2 x lO^^M© cannot make up the entire mass of a 
standard spherical dark halo. Combining these results with those from the 
standard analysis, we find that the halo dark matter may not be comprised of 
objects with masses 2.5 x lO^'^M© < m < 8.1 x 10"^ M0. 



Subject headings: dark matter - gravitational lensing - Stars: low-mass, brown 
dwarfs 



1. Introduction 

If a significant fraction of the dark halo of the Milky Way is made up of Machos 
(MAssive Compact Halo Objects), it should be possible to detect them by searching for 
gravitational microlensing ( [Paczyhski 1986| ; [Petrou 1981| ). As a Macho passes near the line 



of sight to a background star, the star appears to be magnified by a factor 

A = 4^ (1) 

where u = b/vE, b is the distance from the Macho to the line of sight, and the Einstein ring 
radius is given by 

= ^il^^pEi) (2) 

where m is the mass of the Macho, L is the observer-star distance, and x is the ratio of the 
observer-lens and observer-star distances. 

Since Machos are in motion (f ~ = 220km/sec) relative to the line of sight, this 
magnification is time dependent, with A{t) = A{u{t)), where 



u{t) 



"mm ~ 



(3) 



t J 

Here to is the time of peak magnification, and i is the event duration, which can be written 



t = 2rE/v± ~ 130 Jm/M0 days, (4) 
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where v± is the Macho velocity relative to the line-of-sight. For more detailed information, 
see [Paczyhski (1986| ) and priest (1991| ). 



If the halo consisted entirely of objects with masses under about 10~^ Mq the average 
duration of microlensing would be less than 1.5 days, and the events would last only about 
three hours if the halo were made of 10"^ Mq objects. In order to clearly see the shape 
of the microlensing curve for such low mass lenses, images of a lensed star must be taken 
in rapid succession during the event. Such an experiment was undertaken by the EROS 
collaboration ( [Aubourg et al. 1995|) , in which a total of about 82,000 stars were imaged up 



to 46 times per night for several months. No microlensing events were found, and it was 
reported that objects with masses 5 x 10~^ Mq < m < 5 x 10~^ Mq can not comprise the 
entire Galactic Halo at the 90% c.l. We have not followed this approach, but we have used 
a different technique also capable of setting limits on low mass Machos. 

The MACHO collaboration has been monitoring the brightnesses of several million 
stars in the LMC, SMC, and Galactic Bulge since 1992 July using the 50' telescope at Mt. 
Stromlo, Australia. A dichroic beamsplitter and filters are used to provide simultaneous 
measurements in red and blue passbands ( [Alcock et al. 1996b| ). The observing strategy 



for the first two years of LMC data was designed to be sensitive to objects with masses 
m > 10"^ Mq, so a typical LMC field was generally imaged at most once or twice per clear 
night. Therefore microlensing events with durations under a few days will have very few 
magnified points on their light curves and would show up in our data as upward excursions 
of one, two or three consecutive measurements, occurring on stars which otherwise appeared 
completely normal. In this paper we search specifically for such short duration "spike" 
events. Clearly, if any spikes are detected, no conclusion could be drawn regarding their 
origin since there would be insufficient detail in the lightcurves. Therefore, the technique 
described here is most useful when few if any spikes are found, in which case useful upper 
limits can be placed on the prevalence of low mass Machos. After applying selection criteria 
described below, we do not find any such spikes and so are able to strongly constrain the 
existence of low mass objects in the halo of the Milky Way. 



2. Event Selection and Detection 

The analysis reported here uses the first two years of LMC data ( [Alcock et al. 1995a 



Alcock et al. 199"6a| ; [Alcock et al. 1996b|) . Twenty- two fields of 0.5 square degree each were 



monitored on every clear night from 1992 July 20 to 1994 October 26, for a total of 10827 
observations. A total of about 8.5 million stars are used in this analysis. The images are 
taken with a refurbished telescope system (pjart et al. 1996[) and a special purpose camera 
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system ( [Stubbs et al. 1993| ; [Marshall et al. 1994| ), photometrically reduced using a special 



purpose code named Sodophot ( [Bennett et al. 1996| ), and assembled into time series for 
analysis. Each lightcurve consists of many measurements of the flux of a star in two filter 
bands (called "red" and "blue"), as well as estimated errors in the flux measurement and 
several quantities used to detect probable systematic error in the measurement. These 
quantities include the crowding, of PSF fit, missing pixel fraction, cosmic ray flag, and 
sky background. As in the standard analysis these measures are used to remove suspect 
data before any further analysis is performed. 

Also, as in the standard analysis, several properties of the expected microlensing 
signal are used to eliminate stars and events which are unlikely to arise from microlensing. 
After imposing such selection criteria, it is necessary to calculate the number of actual 
microlensing events which would have be removed by these cuts, and this detection efficiency 
calculation is discussed in the following sections. 

Because this search is for short duration events, one of the most powerful signatures of 
microlensing, the shape of the lightcurve (eq's [|I| and 0), cannot be used as a selection 
criterion. Since there are three free parameters for the microlensing lightcurve shape, we 
would need four or more observations during the event to get a meaningful fit. Therefore, 
any phenomenon which causes a significant upward excursion in one or two observations 
could be mistaken for very short duration microlensing. Looking through our data, we 
find many one-observation excursions. A partial list of causes includes satellite tracks and 
glints, telescope slips, and asteroids. In order to reduce this background, we consider only 
those instances where two or three exposures were taken of the same same star on the same 
night. Since each exposure gives both red and blue fiux measurements, we define a "quad" 
as a sequence of 2 (or 3) exposures of a star which are on the same night. In the first two 
years of LMC data we have 1.44 x 10^ quads with two measurements on the same night 
and 5.8 x 10^ quads with three. We then require that all the measurements in the quad 
be significantly magnified, and we also make use of the fact that for very short duration 
microlensing, the measurements of the star on the previous and following nights should not 
have significant upward excursions. Thus each quad represents a potential detection of 
short duration microlensing events, and because the microlensing rate is proportional to 



-1/2 



([Griest 1991[) a substantial number of quads would be magnified if the Milky Way 



Halo were made of low mass Machos. 

We define the magnification of a given measurement as 

A = f/f, (5) 

where / is the fiux of the measurement and / is the median fiux of all the points in the 
relevant pass band that pass the quality cuts described above. In order to reduce the 
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statistical probability of random fluctuations in measured flux giving false triggers, we 
require that all four (or all six) measurements have positive excursions of more than 4a, 
where a is set by the flux measurement error. That is, we require that the magnification 
be above a threshold magnification At, with At — 1 = 4(Tmax, where cxmax is the largest 
magnification error of the four or six quad measurements. By using only these sets of 
measurements and setting a threshold proportional to the error, the probability of false 
event detection can be greatly reduced while still allowing a strong limit to be set. With 
a threshold At — 1 = ^a^sx and a total of 1.5 x 10^ quads, the expected number of false 
triggers from statistical fluctuations is (3.2 x IQ-^)"^ x 1.5 x 10^ = 1.6 x 10"^°. (This 
estimate assumes Gaussian errors, so realistically the non-Gaussian tails will increase this 
number substantially.) The 4crinax threshold was chosen a priori because after investigation 
of several possible analysis methods on a small subset of the data, it appeared most likely 
to give few (if any) events and a signiflcant detection efficiency. 

In addition to the above criteria, we demand that there be a reasonable number of high 
quality data points on the star in order to accurately determine the baseline flux. Thus we 
cut on the number of simultaneous red and blue data points and the average photometric 
error. In our data set we have found a large number of periodic and non-periodic variables 
which can be eliminated due to the fact that they vary continually, so we also demand 
that the lensed star not be flagged as a variable. Next, since microlensing is so rare, we 
do not expect more than one microlensing event to take place on a given star, so we can 
also eliminate stars in which two or more such events are found. Finally, since gravitational 
lensing should magnify all wavebands by the same amount, we can demand that the four or 
six spike points be achromactic within errors. Therefore we make the following set of cuts: 

1) The star must have at least six observations in which both the red and blue data 
points pass the cuts in crowding, seeing, etc. 

2) The star must have V — R< 0.9. 

3) Both the red and blue points in the measurements in the quad, immediately before 
the quad, and immediately after the quad must pass the cuts in crowding, seeing, etc. and 
have a magnification error less than 0.5. 

4) The "robust" (that is, a fit with the highest and lowest ten percent of the 
data excluded) of a fit to a constant flux must be less than 0.9 in both red and blue. 

5) A — 1 > 4(Tinax for all points in the quad, where (Jmax is the maximum error of the 
measurements in the quad. 

6) ^4 — 1 < 4(7inax for the red and blue points in the measurement previous to the quad 
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and the measurement following the quad. 

7) the measurements in the quad must be achromatic within errors, that is A < 2crA, 
where A = |A>Mb - 1|- 

8) There may be at most one event per star. 

Cuts 1, 2, 4, and 8 are for the elimination of false events caused by variable stars, cuts 
1 and 3 are to eliminate background caused by poor photometry, cut 5 is to reduce the 
chance of statistical fluctuations or single-observation glitches causing an event, and cut 6 
is to eliminate longer duration events. See Figure |I| for an example of our data and a Monte 
Carlo event which passes these cuts. 

These cuts were run on the first two years of LMC data, and no events were found. 
To test the robustness of these cuts, the analysis was run with thresholds varying between 
ScTmax and ScTmax- While a few events were found at thresholds below 4(Tmax5 none were 
found at the a priori threshold of 4crmax- This will be discussed further in Section ^. 



3. Theoretical Event Rates 

In order to use our non-detection of spike events to make a statement about the content 
of the Milky Way halo, we need to predict the number of events we would have expected 
to find if the halo is made of low mass Machos. Thus we need to know the efficiency 
with which our experiment, in combination with the above selection criteria, would have 
detected short duration microlensing events if in fact such events occurred. The search for 
spike events is sensitive only to durations of about 0.1 to 4 days, corresponding to masses of 
about 10^^ to 10^'^ Mq (see eq. 0]), and we need to use a halo model to make a connection 
between event duration, rate, and Macho mass. 

We first consider a simple spherical halo model with mass density 

P{r) = Po4i4 (6) 
+ 

where po = 0.008 MqPc"^ is the local dark matter mass density, r is the distance to the 
center of the Galaxy, Ro = 8.5 kpc is the distance of the sun from the Galactic center, and 
a = 5 kpc is the Galactic core radius. [Griest (1991| ) showed that the microlensing rate of a 



5-function mass distribution is given by 

r = 1.60 X IQ-^UT/^Jm/MQ events/year, (7) 

where ut = u{At), and At is the magnification threshold for an event. (The expression 
for r is slightly different from that given in Griest (1991) because we use 50 kpc for the 
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Fig. 1. — A typical Monte Carlo event which passes the cuts used in this analysis. The 
top two curves are an expanded view of the quad in which the event occurred. Note that 
the points in the quad are above the threshold of 1.204, and the previous and following 
measurements are below it. The solid line is the theoretical microlensing curve added to the 
data. On the bottom are the entire year 2 hghtcurves for the same event. The magnified 
quad points are clearly visible at day 540. 
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distance to the LMC, rather than 55 kpc.) In an experiment the number of microlensing 
events which are expected to be detected is given by 

iVe^p = TEE (8) 

where F is the microlensing rate (in event /year/star), E is the effective exposure (in 
star-years) and £ is the average detection efficiency. 

The probabihty of an event occurring with duration imin <i < tmax can be written 

^min 

where, in the approximation of a stationary hne of sight, the distribution of event durations 
is given by ( [Alcock et al. 1996 bj ) 



^ 32utLpo j r%{x)Ae-Q 

di mvii^ J A + Bx + x^ ^ ^ 

^ 

where A = {R^ + a^)/L'^, B = — 2ro cos6cos//L, b and / are the galactic coordinates of 
the source star, vq = 220 km/sec is the solar circular velocity, and Q = Ar'^{x)/{vQp). The 
probability that such an event is actually observed depends strongly on the duration, since 
the spike selection criteria above eliminate the possibility of detecting long duration events 
or events which last only an hour. We can thus define the average efficiency of observing 
such events as 

imax 
^min 

where e{i) is the detection efficiency as a function of event duration. 



4. Finite Source Effects 

When the impact parameter of the lens is comparable to the size of the lensed object, 
the magnification can differ significantly from the point source approximation given in 
equation (|l|). For a lensed star with radius i?*, we define 

U* = ^, (12) 
rEix) 
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as the "effective radius" of the star (the radius of the star normahzed to the Einstein 
radius and scaled to the lens plane). The point source approximation will then break 
down completely for u < because only the fraction of the surface of the star inside the 
Einstein ring radius will be significantly magnified. In the case of a star of constant surface 
brightness, we use ( [Witt fc Mao 1994|) 



for u = U^, and 



2(u - 1 + f/2 



A = j. . , Ui-,n,k 



TT 



u-U, 8 + (M2-f/2) „/7r 



./4 + (m- V2' 



F[i:M (14) 



for u ^U^, where 



n 



An 



V4+(^z-f/,)2' 

and F, E, and 11 are elliptic integrals of the first, second, and third kind. 

There are two effects from the finite source size. First, the maximum possible 
magnification 



_ ;4 + f/2 

^max — 77 (15) 

becomes very low for lower mass lenses. For a star with R = 10 and a lens at x = 0.5, 
we have Amax = 18.3 for a lens with m = 1O~^M0. However, with m = 1O~^M0 we have 
^max = 2.08 and with m = lO'^M© we get Amax = 1.15. Because we are searching for 
significant magnifications, this effect would tend to lower the detection efficiency for lower 
mass lenses. The second effect of a finite source size is that the star is magnified for a longer 
period of time. This occurs because a fraction of the star can be close enough to the lens 
to be significantly magnified even if the lens is far from the center of the star (see Figure 
This effect increases the detection efficiency for low mass lenses whose average event 
duration is shorter than the minimum detectable point source event time scale. These two 
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Fig. 2. — A plot of light curves for various source radii. For larger sources, the maximum 
magnification decreases, but the width of the curve increases. 
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effects mostly cancel each other out, with the detection efficiency increasing slightly for 
very low mass Machos. 

It is clear that the shape of the light curve is strongly dependent on the lens distance x, 
so the detection efficiency is a function of both x and t. Since t is a function of both x and 
v±, we use e = e{x,vj_) in order to simplify the efficiency analysis. The average detection 
efficiency (eq. |]ll|) then becomes 

If f clT 
8 = — I dx I dv±- — ; — ^{x,v±) (16) 
r J J dx dv± 



where Vmax is some upper limit on the perpendicular velocity of the lens. The differential 
rate used in equation ([T^) can be found from equation ( pUj ) using a simple change of 
variables: 

dT _ 2uTLpovli Ae-'^IK 
dxdv± mvn A + Bx + 



5. Monte Carlo and Detection Efficiency 

To measure the detection efficiency e{x,v±), a Monte Carlo simulation was performed 
in which randomly generated microlensing events were added to each star in the database. 
Then the same analysis used to search for spike events was performed on these simulated 
data sets. As a function of the lens position and perpendicular velocity of these events 
one can then find the fraction of simulated events which were recovered and define this 
as e{x,v±). From simple geometry one expects microlensing events to have a uniform 
distribution in minimum impact parameter Mmin- A minimum error of 0.0 14 A is added to 
each data point, so the minimum Sumax threshold is given by At = 1.042, or Wmin = 2.262. 
(As stated previously, the analysis was run on several thresholds varying between Sumax and 
ScTinax, SO the minimum threshold was used to set the upper limit of the Mmin distribution. 
This somewhat lowers the detection efficiency for higher thresholds, but the differential 



microlensing rate given in eq. is correspondingly higher due to the factor of ut and 
there is no net effect in the final result.) Therefore, in performing the Monte Carlo, the 
simulated microlensing events were added with a uniform distribution of Mmin from to 
2.262. To adequately sample the widest possible range of event durations and finite source 
lightcurve shapes, the events were generated using a distribution of x which was uniform 
over < X < 1, and a v± distribution uniform over < f _l < fmax = 667 km/sec. In order 
to improve statistics the simulated events were forced to peak during a quad (that is, to 
between the time of measurements previous to and following a quad). The total exposure 
time used to calculate N^^p was adjusted accordingly by using the total "quad time" rather 
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than the length of the observing run. Thus no simulated events were added during weeks 
when the telescope was down. 

The shape of the light curve is a function of the radius of the star for low mass Machos, 
and the radius of a star is correlated with its magnitude. Because brighter stars tend 
to give lower errors in measured magnification, it follows that the shape of a light curve 
is correlated with the event detection threshold. Therefore the radius of the source was 
estimated from its color and magnitude, and this radius was used in conjunction with the 
observer-lens distance x to determine the shape of the simulated light curves. Although 
limb darkening can change the shape of the light curve ( |Witt fc Mao 199^ ), limb darkening 



coefficients are not well known for such a large sample of stars. An investigation into the 
effect of limb darkening has shown that the resulting uncertainties in magnification are 
much smaller than those caused by the uncertainty in the radius of the source star, so the 
effects of limb darkening are ignored in this analysis. 

Since the detection threshold is proportional to the maximum error of the points in a 
quad, it is important to treat the errors correctly when adding a fake microlensing event to 
a light curve. The error of the magnification A can be approximated by 



a={a', + iO.OUf/fy + f/p)^ (18) 

where as is the error from sky background, / is the total measured fiux of the star, and 
0.014/// is the minimum error added after photometric analysis. The minimum error 
is also added in the standard analysis photometry ([Alcock et al. 1993| ; |Alcock et al. 



1996 b| ; [Alcock et al. 1995b| ) in order to account for several sources of systematic error 



in the photometry. There are two limiting cases of the above formula: "sky-dominated" 
error and "fiux-dominated" error. In the case where the error comes mostly from the sky 
subtraction, the error in the fiux does not change significantly when fiux is added to create 
a simulated microlensing event. On the other hand, when the error is dominated by the 
Poisson statistics of the fiux, the error in the added fiux should be scaled by ^/f . Because 
the fiux error is larger, giving higher thresholds, this case was used in the Monte Carlo to 
get a conservative measure of the efficiency. (The minimum error was subtracted before 
the errors were adjusted and then added in again afterward, just as the minimum error is 
added in after the photometric reductions in the standard analysis.) The true efficiency is 
expected to be closer to the fiux error case, because in general the sky error contributes 
significantly only to dim stars which already have a fairly low efficiency due to their large 
error bars (large threshold). 

Because of our crowded fields and poor seeing, many of our photometric objects 
are actually blends of two or more stars. When a blended star is lensed, the measured 
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magnification can be significantly smaller than the true magnification. To quantify 
this effect, artificial stars were added to real images taken under a variety of observing 
conditions and the photometry code was run again. We thus created a series of response 
functions of recovered vs. added flux, and when a simulated microlensing event is generated 
the photometric object is matched to one of these response functions using the object's 
magnitude. The observing conditions on each point in the light curve are then matched to 
similar conditions in the response function, and the recovered flux is added to the data. A 
detailed description of this analysis can be found in [Alcock et al. (1996b| ). 

The effects of blending on the detection efficiency are twofold. First, the lower 
measured magnification will lower the efficiency as fewer quads will be above the threshold 
magnification. However, if an object is a blend of two or more stars, then there are two 
or more stars that may be lensed. Thus the total number of stars in the data set is larger 
than the number of photometric objects, and our total exposure is significantly increased. 
A Monte Carlo event with a large blend fraction (and significant finite source effects) can 
be seen in Figure ^. 



6. Results 

Because the blending effects described above are a function of the magnitude of the 
lensed object, it follows that the detection efficiency depends on the magnitude of the source 
as well as the lens position and velocity. Using an infinitesmal bin width, the efficiency 
e(x, f _L, M) can be written 

dN^ec /dxdv±dM ^ ^ ^ ^ 

dNs^dd/ dxdv±dM 

where M is the average of the red and blue magnitudes of the lensed star (this value is 
the unblended stellar magnitude which is determined by the blending response function), 
A^add is the number of fake events added to the data, and Aj-ec is the number of these events 
recovered by the analysis. To calculate the number of expected events, one must integrate 
the efficiency over the stellar luminosity function n{M): 

A'exp = J dM J dx J dv±-^^^^e{x,v±,M)n{M)T (20) 



where T is the effective exposure, or "quad time" of a given star. To calculate the efficiency, 
we have 

dN 

^ _ ^^)S^y^ _ ^^,)5(M - M,) (21) 



dxdv±dM 
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Fig. 3. — A Monte Carlo event with significant blending and finite source effects. The shape 
of the theoretical light curve differs from a simple point source light curve, and the magnified 
points are below the added magnification because only a fraction of the object is lensed. The 
event still passes the cuts, and the magnified quad can be seen at day 450 on the entire year 
2 light curve. 
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where the sum is over all recovered events. Because we add one simulated event to each 
object in the database, we can write 

^ = nUM) (22) 

where nsod(^) is the Sodophot object luminosity function (that is, the luminosity 
distribution of objects recovered by the photometric reductions of the images). This 
function is used because the response function stellar magnitudes follow the Sodophot 
object distribution, and the response functions are chosen uniformly from this distribution 
( lAlcock et al. 1996b| ). The events are added uniformly in x and v±, which gives 



-HsodiM) 



dxdv±dM dxdvj_ 

= const. X ?isod(M) 
risodiM) 



(23) 



where the factor 1/fmax is a normalization constant. The number of expected events then 
becomes 



dxdv\ 



nsod(M, 



-^i^max (24) 



Approximately 6.5% of the stars used in this analysis can be found in more than one 
field (Alcock et al. 1996b), and we must scale the number of expected events accordingly. 
However, the possibility of double counting in this analysis occurs only when double or 
triple exposures are taken on overlapping fields on the same night. This only happens in 
about 2/3 of our quads, so we subtract 4.4% from our number of expected events rather 
than 6.5%. 

A plot of A'^exp vs mass for a (5-function mass distribution is given in Figure ^. At the 
peak at 1O~^M0, about 17 events would be expected to have been found if the halo is as 
modeled in equation (P), and consisted entirely of Machos of that mass. Equivalently, we 
may convert the number of expected events into upper limits on the allowed halo mass 
fraction that can be contributed from objects in the excluded mass range. Such a plot is 
given in Figure ^ Using the fact that no events were found we can place strong limits on 
the halo of the Milky Way. The 95% Poisson c.l. when A^^obs = is A'exp = 3 events, so for 
the simple spherical halo model masses between 2.5 x 10^'^ M© and 5.2 x 1O~^M0 are ruled 
out at the 95% c.l. Although these limits are for a 5-function mass distribution, any model 
distribution containing a combination of masses in this range is also ruled out at the 95% 
c.l. ( IGriest 1991]; [Alcock et al. 1996b|) . 
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Macho Mass (M^) 

Fig. 4. — A plot of the number of expected events vs mass for a 5-function mass distribution. 
With no events found the 95 % c.l. upper hmit is 3 events, and the region of the curve above 
this hmit are excluded. Also shown is the number of expected events from the EROS CCD 
experiment. 
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Macho Mass (Mq) 

Fig. 5. — A plot of allowed halo mass fraction vs mass for a 5-function mass distribution. 
The region above the solid hne is excluded at the 95% c.l. Also shown is the halo fraction 
upper limit from the EROS CCD experiment. 
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Also shown in Figures |^ and ^ are the results from the EROS CCD experiment 
( |Aubourg et al. 1995| ). Although they give a stronger limit for m < IO^^Mq, the limits set 
by this analysis give the strongest limits to date for 1O~^M0 < m < lO^^M©. 

As mentioned previously, the analysis was run with thresholds varying from 3o"max 
to 5o"max in order to determine the robustness of the analysis. A plot of the number of 
expected and observed events as a function of threshold can be found in Figure ^, and 
the corresponding plot of halo fraction upper limit vs. threshold can be found in Figure 
1^. A total of 11 events were found which passed various thresholds between ScXmax and 
3.75(Jmax, but no events were found at thresholds of 4(Tmax and higher. Of the 11 events 
found, 8 were on stars with V < 17.5 and are likely low level variables which fell through 
the variable star cuts. Of the three remaining events, two occurred in the same field on the 
same night which indicates possible problems with the observations, and inspection of the 
second image in the quad shows a likely telescope slip during the exposure. The remaining 
event passes only the Samax threshold cut. In order to reduce these backgrounds in future 
analysis runs, stars with V < 17.5 will be cut as will images with more than one event, and 
it has been determined that these cuts will reduce the number of expected events by about 
20%. However, because neither the number of expected events or the number of observed 
events varies drastically with threshold, the analysis using the a priori 4(Tmax threshold is 
robust and the limits set on the abundance of low mass Machos are valid. 



7. Combined Analyses 

The standard analysis method of fitting microlensing curves to the data is sensitive to 
Machos of masses 10~^ M© < m < 1 Mq, and it would be useful to combine the results of 
the two types of analyses. To avoid double counting of events which could pass both the 
spike and standard cuts, we ran the standard analysis cuts on any simulated events passing 
the spike event cuts, and any events passing both sets of cuts are so flagged. The efficiency 
is then recalculated with the flagged events considered as failing the cuts, and the number 
of expected events can then be added to the number of expected events from the standard 
analysis. However, when adding the number of expected events the number of observed 
events must also be added. The standard analysis of the first two years of LMC data yields 
eight likely microlensing events ( [Alcock et al. 1996a| ) which would make the combined limit 
on halo fraction very weak. However, the eight events all have durations t > 34 days, and it 
is very unlikely that these events were caused by machos with m < 0.1 Mq. Therefore the 
standard analysis efficiency was recalculated with a cut such that the event duration must 
be shorter than 20 days. (20 days was chosen in order to be conservative.) The number 
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Fig. 6. — Number of expected events as a function of threshold for three values of Macho 
mass. Also shown are the number of observed events as a function of threshold. 
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Fig. 7. — Halo fraction upper limit (95% c.l.) as a function of threshold for three values of 
lens mass. The upper limit rises at lower thresholds because events are detected at these 
thresholds and the 95% c.l. upper limit on A^exp increases accordingly. 
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of expected events thus drops significantly at m > O.IMq, but we also have no observed 
events and are able to place strong limits on lower mass objects. The number of expected 
events and halo fraction vs lens mass can be found in Figures |^ and ^. Here it can be seen 
that Machos of masses 2.5 x 1O~^M0 < m < 8.1 x 1O~^M0 can not make up the entire 
halo mass, and lenses in the range 1.88 x 10^^ M© < m < 2.5 x 1O~^M0 comprise at most 
20% of the halo dark matter. 



8. Power Law Halo Models 



The halo models of Evans ( [Evans 1994| ; [Evans fc Jijina 1994] ) allow for rising or falling 



rotation curves, fiattened halos, and various disk contributions to the total galactic mass. 
These models are also called "power law" models because at large galactic radii R, the 
circular velocity fcirc for some model parameter /3. The parameters used to describe 

the mass and velocity distributions in these models are as follows: 

P At large galactic radii, (3 = gives a flat rotation curve, /3 < gives rising curve, and 
P > gives falling curve. 

q Halo flattening parameter, q = 1 gives spherical halo, g = 0.7 represents ellipticity of 
E6. 

Vo Normalization velocity. 

Rc Galactic core radius. A large Rc gives a massive disk. 
Rq Radius of solar orbit. 

The differential event rate dV /dxdv^ can be derived as a function of these parameters 
( [Alcock et al. 1995c[) , and it is then straightforward to calculate the number of expected 
events using equation ([2^) . Limits were calculated for the same models used in [Alcock et al.[ 
[(1996b[) , and the parameters used are found in Table Also shown for each model in Table 
|I] is the total mass inside 50 kpc from the center of the Milky Way, which we call M50. 

Model S is the simple standard spherical halo described in Section ^, and model A is 
the power law model equivalent. Model B has a rising rotation curve and a more massive 
halo, while model C has a falling curve and a less massive halo. Model D has a flattened 
halo, and models E, F, and G have more massive disks. Model E has an extremely massive 
disk and a very light halo, and this model is probably inconsistent with estimates of the 
mass of the Milky Way. 
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Macho Mass (M^) 

Fig. 8. — Number of expected events as a function of Macho mass after combining the 
standard and spike analyses, with a cut on events with t > 20 days. The spike result is 
plotted with the dotted line (with the double counted events subtracted), the standard 
result with the dashed line, and the combined result is shown with the solid line. 
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Fig. 9. — Halo fraction upper limit (95% c.l.) for the combined spike and standard analyses 
and with t < 20 days. 
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The number of expected events was calculated for these models and then combined 
with the results from the standard analysis as described in Section |^. In Figure |1^ we plot 
the resulting number of expected events as a function of lens mass for a (5-function mass 
distribution using the simple spherical halo model and the seven power law halo models 
described above. Figure |ll| is a plot of allowed halo fraction vs mass for the same models. 
For the models with more massive halos, only about 30% of the halo can be comprised of 
Machos in the range of 9.5 x 10"'' < m < 2.9 x 1O~^M0. The limits get weaker for less 
massive halos, and little useful parameter space can be excluded for the extreme "maximal 
disk" model E. 

The differences between the limits among the various models is primarily because the 
number of expected events is directly proportional to M50, or the number of Machos in the 
halo. We can get more model-independent limits by removing this factor and plotting the 
total allowed halo mass from Machos inside 50 kpc, rather than the halo mass fraction, 
as a function of mass. This plot is shown in Figure |12[ Ignoring the unlikely model E 
we see that, independent of model, no more than 10^^ Mq of the halo mass inside 50 kpc 
can come from objects of mass 1.85 x 1O~^M0 < m < 6.5 x 10~^Mq, and objects in the 
range 3.2 x 10"'' Mq < m < 1.87 x 10^^ Mq can not make up the entire canonical value of 

4.1 X 10" Mq. 



9. Conclusion 

We have extended the sensitivity of the MACHO experiment to two orders of magnitude 
lower in mass using existing data and without changing observing strategy. Objects with 
masses 2.5 x 10"'' M© < m < 8.1 x 10~^Mq (roughly one Mars mass to 80 Jupiter 
masses) can not comprise the entire standard spherical halo mass, and Machos in the range 
1.88 X IO^'^Mq < m < 2.5 X 10"^ Mq make up less than 20% of the halo. Independent of 
halo model, objects in the range of 3.2 x 10^'' Mq < m < 1.87 x 10^^ Mq can not make up 
the canonical halo mass inside 50 kpc of 4.1 x 10^^ Mq, and less than lO^'^ Mq of the halo is 
made from Machos with masses 1.85 x 10~®Mq < m < 6.5 x 10"^ Mq. These limits are the 
strongest published to date. 
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Fig. 10. — A plot of the number of expected events vs mass for the standard model and 
seven power law halo models. The results shown are for the combined spike and standard 
analyses, and with i < 20 days. The line at N^xp = 3 is the 95% c.l. upper limit, and the 
regions of the curves above this line are ruled out. 
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Mass (Mq) 

Fig. 11. — Upper limits on Macho fraction of the halo vs lens mass for the spherical and 
power law halo models (line coding is the same as in Figure 10). The results shown are for 
the combined spike and standard analyses, and with i < 20 days. The regions above the 
curves are ruled out at the 95% c.l. 
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Fig. 12. — Upper limits on the total mass of Machos interior to 50 kpc as a function of 
lens mass. The results shown are for the combined spike and standard analyses, and with 
i < 20 days. The regions above the curves are ruled out at the 95 % c.l. Objects of mass 
3.2 X 10"''' Mq <m< 1.87 x 10"^ M© can not make up the canonical value of 4.1 x 10^^ M©, 
independent of the model used. 
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Table 1. Power Law Halo Model Parameters 



Model 




Q 


vo (km/sec) 


Rc (kpc) 


Ro (kpc) 


Mso (IO^^Mq) 


S 








5 


8.5 


4.13 


A 





1 


200 


5 


8.5 


4.62 


B 


-0.2 


1 


200 


5 


8.5 


7.34 


C 


0.2 


1 


180 


5 


8.5 


2.36 


D 





0.71 


200 


5 


8.5 


3.74 


E 





1 


90 


20 


7.0 


0.82 


F 





1 


150 


25 


7.9 


2.10 


G 





1 


180 


20 


7.9 


3.26 



